A domain wall model for spectral reflectance of plant leaves 
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We model a plant leaf by using two-dimensional domain walls with internal structures. Such domain 
walls can be found as soliton solutions in field theory describing magnetic materials. The radiation 
scattered by such domain walls behaves quite similar to the spectral reflectance of plant leaves. The 
model nicely simulates the spectral reflectance of a plant leaf as a function of the wavelength. 
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/. Introduction. The study of the spectral behavior 
of vegetables is often used to investigate the character- 
istics of the electromagnetic radiation reflected by plant 
leaves, a plant as whole or a collection of plants. We 
should emphasize that the study of the spectral behavior 
in general means the study of spectral reflectance, spec- 
tral transmittance and spectral absorbance. Concerning 
the interaction of the solar energy with vegetables, the 
plant leaf plays a special role because it basically realizes 
the photosynthesis that is responsible by the formation of 
carbon compounds. The whole configuration of a plant 
leaf including form, position, structures, etc, adapts to 
receive in a most efficient way the solar beams, air and 
water which are necessary for realizing photosynthesis. 
Thus the knowledge of the properties of a plant leaf is 
clearly fundamental, especially in the study of the re- 
flectance of a plant or even a culture of plants. Such 
predominance of the plant leaf is so important that the 
area of the other part of the plant in contact with the 
solar radiation is completely neglected Q, Q • It is obvi- 
ous, however, that the data obtained from a unique plant 
leaf cannot be used directly, without modifications, to a 
culture of plants. This is because there exist qualita- 
tive and quantitative differences between the aspects of 
unique plant leaf and a culture of plants in a field. The 
structures of the cells that constitute the three tissues 
of the leaves vary with the specie and ambient condi- 
tions. The compounds of the leaf that are of consider- 
able importance in the study of interaction of a leaf with 
the radiation arc: cellulose that are found in the cellular 
walls, solutes like ions and molecules, intercellular spaces 
and pigment inside chloroplasts such as carotene, xantho- 
phyll and chlorophyll. The spectral behavior of a plant 
leaf is a function of its compounds, morphology and in- 
ternal structure. Since the characteristics of the plant 
leaves are genetically controlled there exist differences in 
the spectral behavior among groups that are genetically 
distinct. The Fig. [I] shows possible paths of the radia- 
tion or incident energy over a plant leaf. A small amount 
of radiation is reflected from the cells around the sur- 
face of the plant leaf; the major part of the radiation is 
transmitted to the spongy mesophyll, where the cellular 
walls for incident angles sufficiently large reflect the in- 



cident radiation rays These multiple reflections are 
essentially a random process where the rays change their 
directions inside the plant leaf. Because the large number 
of cellular walls, some rays go back toward the source of 
radiation, while others are transmitted trough the plant 
leaf. The thickness of the plant leaf clearly affects the 
transmittance of the radiation, being larger than the re- 
flectance for thin leaves. The opposite happens for thick 
leaves. The spectral response of a target, as a plant leaf 
considered in the present work, deal with the graphic 
representation of the reflectance in narrow and adjacent 
bands of wavelengths. This representation gives us the 
effect of the interaction of incident radiation and the tar- 
get under investigation. Thus, the amplitude variations 
in spectral response give us clues about the spectral prop- 
erties of the objects [3,|5|)l10- The curve of reflectance 
of a plant leaf can be divided into ultraviolet region, visi- 
ble region (400 - 700 nm), near infrared (700 - 1300 ran), 
and mid infrared (1300 - 2600 nm). However, the ultra- 
violet region usually is not considered because the most 
part of this radiation is scattered away and the vegetables 
does not make use of this radiation. The mid infrared is 
not considered here because of limitations on the sen- 
sor equipment in this region. The curve of reflectance 
obtained experimentally is shown in Fig. 2. 

In this work, we model the plant lea f by a domain wall 
with internal structure 0, & E3, E3, E3 In the follow- 
ing section we show with details how this can be possible. 
The model is mainly based on the reflection probability 
that light neutral particles can suffer as they collide with 
such domain walls — in Ref. Q this issue was initiated 
in another context. In our model the light neutral par- 
ticles play the role of the electromagnetic radiation and 
the domain wall with internal structure plays the role of 
a plant leaf. 

The internal structures are small domain walls that 
we identify with the spaces among the spongy meso- 
phyll. Domain walls are soliton solutions that appear 
in systems that exhibit spontaneous symmetry breaking, 
such as the Higgs mechanism applied in Particle Physics 
and the Ginzburg-Landau theory applied in Supercon- 
ductivity. These systems usually present non-linearity 
into equations of motion. Such non-linearity is respon- 
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such that V(4>,x) can be written as 




FIG. 1: Transversal section of a plant leaf with possible paths 
of radiation rays. (1) upper epidermis, (2) palisade mesophyll, 
(3) spongy mesophyll and (4) lower epidermis. 



sible for the formation of soliton solutions whose energy 
density is localized around a region of the space, and 
remains stable all times. 
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FIG. 2: The spectral reflectance of a peanut leaf in the visible 
and infrared regions obtained experimentally. 

77. The domain wall model. The dynamics of a system 
with two coupled real scalar fields <p an d X is described 
by a field theory with the Lagrangian 

c = hc> ll ci>d»tt> + hd^x-v{<t>,x)- (i) 

Where the potential V(4>, y) belongs to a class of soliton 
models easily integrable 0, U 03 ED, 13 

given by 



V 



1 fdW\ 2 1 fdW\ 2 



2 V 



2 V d<t> J ' 



(2) 



where W is known as "superpotential" that for our pur- 
poses here it is sufficient to make the following choice 



H ' = "1 y-a 2 0) +^X 2 



1 

— ! 

2 

1 2 4 

2"- 



V{4>, X ) = ^ 2 W> 2 -a 2 ) 2 + (2 M 2 + ^)</>V -vna\ 2 



(4) 



The scalar field <j> describes the domain wall that we iden- 
tify with a plant leaf. The scalar field y describes the 
internal structures that appear as domain walls, which 
can form disk like structures |13|| . These domain walls 
will be identified later with cellular walls inside the plant 
leaf. The meaning of the parameters is, a and fi will be- 
come clear later. The field equation of motion are given 
by 



and 
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Now we make use of perturbation theory, up to first or- 
der, around the soliton solutions <f> and y into the equa- 
tions of motion (JSJ and ©, as x — X + £ an d <jj = 4>, 
where £ represents the relevant fluctuations for the inter- 
nal structures. Considering such fluctuations having the 
form £ = £(z)e~ 1 ( LOt ~ kj;X ~ kyV ' > the relevant equation for 
the fluctuations reads 



+ k 2 x + k 2 y \az) + 0£(z) = 0, (7) 



which is a Schroedinger like equation. The bar means the 
function is evaluated at the original soliton solutions. By 
considering a domain wall approximately plane and static 
it suffices to work only with a spatial dimension such 
that the fields are functions as <j) = <p(z) and x = x( z )- 
The equations of motion have the following well-known 
solutions [HQ3 



4> = — atanh(^az), x = 0, 



(8) 



6 = — a tanh (2/iaz), x = ^ sza \ / 2 sech (2//az).(9) 

V m 

Where z is a coordinate that is transversal to the domain 
wall. The type I solution JSJ represents a domain wall 
without internal structure while the type II solution @ 
represents a domain wall with internal structure. Both 
solutions have the same energy (or rest mass) Mw = 
aw Aw with surface density 



(10) 



that can be found by taking the difference of the su- 
perpotential evaluated at the vacuum solutions </> = 
±a and x = , i.e., a w = \AW\ - see references 
In the type II solution, as the 



,„< ±a and y = , i.e 

(3) h i mm m m- 
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field (j> — > in the interior of the domain wall, the field 
X develops a nonzero maximum value that characterizes 
its localization into the domain wall. In such regime the 
potential V(4>,x) given in approaches effectively the 
potential 



eff 



¥ (* 2 - 7 « 

2 V » 



(11) 



The effective Lagrangian that now governs the dynamics 
inside the domain wall is given by 



C eff = ^xd^X-Veffix)- 



(12) 



Given the potential Qllfl. this theory is able to describe 
internal structures such as other domain walls given by 
the following solution 



X = ±a 



tanh 



(13) 



where Xi is some tangential coordinate along the domain 
wall that is going to be identified with a plant leaf. The 
structures have energy (or rest mass) Ms = as As with 
surface density given by 



Here we used the fact V — A S, being S t he domain 
wall thickness. Thus 0= (n/2)y/l7 - 4A 4 /A^, and since 
k z = 2ir/\ and S — 1/2/ia we find 7 = 2ir 2 5/A. Where we 
have approximated the area of the internal structures as 
the area of a disk like structure As 27rL 2 and consid- 
ered that the typical size of the structures is related to 
a critical wavelength as L Ao . On the other hand, by 
considering Aw w 2irN\ 2 we define Ayy/N w 27rA 2 as 
the smallest area of radius A covered by a radiation with 
wavelength A. 
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This can be found by constructing the "effective super- 
potential" W e ff — M (is — /7° 2 ^) t na ^ produces the 
effective potential given by (|llf> . As in the previous 
case, the surface density is found by evaluating the su- 
perpotential at vacuum solutions x — ± a \ —■> sucn that 

y A* 

a s = I A W e // 1 . 

Now using the type II solution into equation J7J) we 
find 



(z) + k 2 z t{z) ~ U{z)i{z) = 0, 



(15) 



where — k 2 = —uj 
of the particle momentum 
mass of the particles far from the domain wall, and the 
Schroedinger potential is given by 



m 2 is the z-component 
= 4/i 2 a 2 is the squared 



U(z) 
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(16) 



The reflection probability related to the potential U(z) 
is well-known [l4| 
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with (3 = (tt/2)VT7 -4*V/i and 7 = nk z /2fia. How- 
ever, by noting that the rest mass of the domain wall 
can be related to the rest mass of the internal structures 
as Mw — NM S , we can write 



<JwVw = NasVs — > A 



w - 
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FIG. 3: The spectral reflectance of a plant leaf according to 
the domain wall model, for Ao = 400 nm and 5 — 30 nm. 

III. Results and discussions. The Figures [3 and 
show the reflectance curves obtained experimentally via 
spectrum-radiometer LI-1800 from LI-COR and the- 
oretically via domain wall model, respectively. As we can 
see in these figures, the model (Fig. |3| nicely simulates 
the real spectral behavior of a plant leaf (Fig. [2J. In 
the visible region (400 a 700 nm) one observes two valley 
around 480 nm and 680 nm that evidences the absorp- 
tion of the chlorophylls "b" and "a" , respectively, and 
a peak around 555 nm that characterizes the color of a 
green plant leaf. In the near infrared region (700 nm - 
1300 nm) the reflectance is almost constant. Let us now 
compare the results obtained experimentally with the re- 
sults obtained theoretically via domain wall model, by 
entering with some known parameters of the model. We 
can use the formula below to estimate the critical wave- 
length 

(19) 

For structures (or spaces) with typical dimensions about 
4-12 \im of width and 10-14 /im of length, as the regions 
separating spongy mesophylls, we find critical wavelength 
Ao ~ 400 nm. The parameter S is associated with the do- 
main wall thickness, and can be regarded as an effective 
thickness of the plant leaf, where effectively reflections 
may occur, i.e., the interior region of the plant leaf. The 
Fig. |3 shows the theoretical behavior of the spectral re- 
flectance of a plant leaf according to our model. Un- 
like the spectral reflectance obtained experimentally (see 
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Fig. |5J) the minimal reflectance is zero because only in- 
ternal effects are considered in the model. For the sake 
of simplicity, we disregard reflections from the surface of 
the plant leaf. Note that the spectral reflectance from 
internal structures of plant leaves is well approximated 
by our domain wall model. Although the second point 
of absorption around 565 nm is a bit far from observed 
value 680 nm, the first point of absorption at 480 nm is 
in perfect agreement and the reflectance peaked around 
540 nm is pretty close to the observed value around 555 
nm that characterizes the color of a green plant leaf. 

IV. Conclusions. The model nicely simulates the spec- 
tral reflectance of a plant leaf as a function of the wave- 
length, i.e., its spectral behavior. The domain wall model 
works well for structures in the scale of micrometers {um 
) . This is also the scale of spaces between cells ( « cellu- 
lar walls thickness) , that we can identify with the domain 
wall like structures in our model discussed above. In or- 
der to simulate better the spectral behavior of a plant 



leaf we expect to improve the model in a separate work. 
Such improvements may be concerned with the number 
and the type of domain walls like structures formed inside 
the domain wall that simulate the plant leaf. This may 
improve the description of several intercellular spaces giv- 
ing the real complexity of the internal structure inside 
a plant leaf. Finally we should emphasize that our in- 
vestigations concerning a domain wall model to simulate 
spectral reflectance here applied to a plant leaf can eas- 
ily be extended to a class of planar systems with internal 
structures. 
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